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Recent BaBaR data on the pion transition form factor, whose dependence is much steeper 
then predicted by asymptotic Quantum Chromodynamics (QCD), have caused a renewed interest in 
its theoretical description. We present here a formalism based on a model independent low energy 
description and a high energy description based on QCD, which match at a scale Qo- The high 
energy description incorporates a fiat pion distribution amphtude, 4>{^) = 1, at the matching scale 
Qo and QCD evolution from Qo to Q > Qo- The flat pion distribution is connected, through soft 
pion theorems and chiral symmetry, to the pion valance parton distribution at the same low scale 
Qq. The procedure leads to a good description of the data, and incorporating additional twist three 
effects, to an excellent description of the data. 
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I. INTRODUCTION 



The implications of recent data by the BaBar Collaboration [1] on the transition form factor 7*7 — tt", in our 
understanding of the structure of the pion are being widely discussed These results have cast doubts on 

I '' our understanding on the behavior, as a function of the light-cone momentum fraction x, of the pion distribution 
amplitude (ttZ? A) (ijj. (a;) 0,01 ■ Some investigations have produced scenarios where the pion ttDA is flat, i.e. a constant 
value for all a: |3[, and are in good agreement with the data for the form factor. These scenarios arc compatible 
with QCD sum rules and lattice QCD d, ^ calculations which provide values for the second moment of the 
t:DA which are large compared with the asymptotic value 6x{l — x) Q ^- Model calculations, Nambu-Jona-Lasinio 
(NJL) jllrtlSij and the "spectral" quark model [l^ give a constant pion distribution amplitude, 4'{x) = 1. The pion 
transition form factor calculated in these models, however, overshoots the data [l6j . 

In here we present a formalism to calculate the pion transition form factor (ttTFF) which provides an excellent 
, description in the whole range of experimental data. The formalism consists of three ingredients: i)a low energy 
' description of the ttTFF; ii) a high energy description of the ttTFF; Hi) a matching condition between the two 
descriptions at a scale Qo characterizing the separation between the two regimes. For the low energy description we 
■ take a parametrization of the low energy data to avoid model dependence at Qq. The high energy description of the 
' ttTFF, defined by the pion Distribution Amplitude (nDA), contains Quantum Chromodynamic (QCD) evolution 
from Qo to any higher Q, the recently introduced mass cut-off procedure [2], and ultimately additional higher twist 
effects. The relation between the leading twist (ttDA) and the pion Parton Distribution [ttPD) allows us to use a 
scheme for the determination of Qo that has proven successful to describe parton distributions of mesons and baryons 
' [T7l - [23| . The scheme consists in evolving the second moments of the parton distributions from a scale where they are 
M \ experimentally known to a scale Qo where they coincide with those determined by a low energy description. Once 
Qo is known, higher Q results for any other observable are obtained by evolving the low energy description of the 
corresponding observable from Qo to the required scale Q. Our analysis will systematically explore the region of 
validity of the scheme and find implications of new physics in the intermediate energy region. 

In section II we define the ttTFF both for low and high energies and establish the matching condition at the 
hadronic scale Qo . The matching condition of section II allows to determine the mass parameter of the formalism 
once the hadronic scale is known. Section III is devoted to establish a relation between the ttDA and the ttPD. 
This relation allows to characterize a unique matching scale Qo. Section IV is dedicated to analyze the stability of 
the parameters and the description of the data. We find reasonable agreement with the experimental data for Q^ 
below 15 GeV^ with parameters whose values are to be expected from physical arguments. Two matching scales are 
investigated in the analysis. We realize that none of them is able to reproduce accurately the slope of the data in 
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the intermediate region. In section V we incorporate in a phenomenological way additional twist three effects. They 
lead to a modification of the matching condition. Their overall effect is small at high Q but they are important in 
reproducing the slope in the intermediate region. In this way we are able to get an excellent agreement with the data. 
In section VI we draw some conclusions. Finally, in the Appendix we present a model calculation which overshoots 
the data and therefore signals the importance of evolution in the description of the experiment. Morevoer, it also 
serves to understand the order of magnitude of some of the parameters used in our approach. 



II. THE PION TRANSITION FORM FACTOR. 



The form factor F^t^-Tr relating two photons (real or virtual) to the pion is a very important object of study 
in exclusive processes in QCD. In particular it has been used by comparing perturbative QCD predictions with 
experimental data to obtain information about the shape of the 7rDA0, HI, Experimentally for small virtualit" 
of one of the photons it has been measured by the CELLO [l^, by the CLEO 27 1 and recently by the BaBaR [l| 
collaboration. The last results are in disagreement with previous theoretical expectations. 

In lowest order of perturbative QCD, the transition form factor for the process tt*^ — 77* is given by 



S Jo X 



(1) 



where = —q^, is the momentum of the virtual photon, (j>Tr (a;, Q^) is ttDA at the scale and /^r = 0.131 CeV. 
In this expression, the dependence appears through the QCD evolution of the ttDA. The ttDA can be expressed 
in terms of the Gegenbauer polynomials [2^, [2^ , 

Mx,Q^)=xil-x) anCl/\2x-l)(log-^] , (2) 



n{even)—0 



A? 



where 7„ are the anomalous dimensions 




-■• = ^l'-^''>.^-(;m)(^T2yj- 

/3 — ^^^"^ — -y^ is the beta function to lowest order and Cp — ^2Nc ' ^^^^niii^S that the a„ coefficients are known, 
Eq. ([1]) gives the ttTFF for any value of Q high enough compared with Aqcd ■ The a„ coefficients are obtained if we 
know the wDA, (fyj, [x, Qo)' ^t some momentum Qo, using the orthogonality relations of the Gegenbauer polynomials 

In order to prevent divergences in the integrand of Eq. ([1]) it has been assumed that (pj^ {x) vanishes for x — Q. There 
are no fundamental reasons for this restriction, moreover, several chiral models predict for massless pious ipTrix) = 1 
We use this model prediction as input at Qq in Eq Q and we obtain for the expansion of the ttDA 



This expression has the right asymptotic behavior, limQ_i.oo 4>-n{x^Q) = 6x(l — x). Nevertheless, despite its appear- 
ance, the T^DA defined in Eq. ([5]) has the properties limj.^o+ 07r(a;, Qo) = hm^-^i- 4)t^(x^Qq) ~ 1. In order to cure 
the divergences in the integrand of Eq. ([T]) we follow the proposal of Radyushkin in section II. B of ref. 0, and we 
introduce a cutoff mass M in the above definition of the ttTFF, in terms of the ttDA, 

Q'FiQ') = f -^M^, Q')- (6) 

-J Jo x+ YW 
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We will discuss later the determination of the M. 
At low we can parametrize the data as (30l |. 



F{Q^)^F{Q) 



Q 

a 



■tt'-' 



(7) 



with F (0) — 0.273(10) GeV^^ and a = 0.032 (4) , constants which have been obtained from the experimental study 
of 7r° — 7 6"'" e~. The ttTFF at these low values of can be described by several hadron models. For instance, the 
parametrization can be well understood studying the triangle diagram (see Fig. [T}. This calculation is straightforward 
(see appendix) leading, in the case of massless pious, to F{0) = \/2/ {^i-K^f-^) — 0.273 GeV^^, which coincides with 
the experimental value. The parameter a can be linked to the constituent quark mass, obtaining ruq = m^j \JV1 a ~ 
0.220 GeV. 

Despite the fact that several models reproduce the low energy data, in order to have a model independent expression 
for the form factor, we adopt a monopole parametrization of the -kTFF in the low energy region: 



We have thus introduced two descriptions for the ttTFF, one in the low energy regime defined by F'"^{Q^), and one 
in the high energy regime, defined by Eqs.® and ([S]). Both descriptions will be matched at Qqi the scale previously 
introduced which we now define, li Q > Qo the pion form factor should be calculated using Eqs. ([5]) and (|6]) , and 
for Q < Qo by Eq. ([5]). For Q ^ Qo both expressions should lead to the same result. In this way Qo is defined as the 
matching point between the low and high energy descriptions. Consistency implies that this matching point must be 
the same for any physical quantity. 

Note that our formalism contains two unknowns, Qo and M. The scale Qo is an important ingredient of our 
calculation p^| - |20| and determines the range of validity of the low energy description used in the calculation, i.e. for 
Q < Qo 07r(a;) = 1 and F^^ (Q^) is assumed to be a good representation of the theory. Qo is fixed by the value of 
^QCD, and the physical results should not show any dependence on it. M contains not only the effects associated 
with the constituent quark mass but, as pointed out by Radyushkin also the contribution of some mean transverse 
momentum. 

We next equate the pion form factor calculated via equation ([S]) at Q^ = Qq, where 07r(x, Qo) = 1, 



Qo^(Qo) - ^ / — ^ - ^'"^w^' 



to the value given by the monopole parametrization Eq. ([8]). 
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This is the first equation relating the two parameters. We proceed in the next section to find a second relation 
relation which allows a full determinations of the parameters. 
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III. THE PION DISTRIBUTION AMPLITUDE AND THE PION PARTON DISTRIBUTION. 

In this section we establish a relation between the ttDA and the pion valence parton distribution [ttVPD) which 
will be instrumental in the determination of Qq. 
The ttDA is defined by, 

<5,,zV2/, = J ^e^^- (0|V;(0)7+75TV(An)|^^"(P)). (11) 

In order to introduce a convenient expression for the valence parton distribution let us define the quark distribution 

HI'' i^) = lj ^e^'" iP)\^ (0) 7+ i (1 + ar^) (An) \7r^ (P)) = 1 {6,,, H'^" (x) + ie,,,aH'-^ {x)) (12) 

with a = ±1 for q — u,d and H^^^'^ [x) represent the isospin decompositions. The quark distribution is only non- 
vanishing in the interval — 1 < a; < 1. This quark distribution is related with the standard q (x) and q (x) distributions, 
by 

H^^ ^q{x)e{x)-q{-x)e(-x) , (13) 

with 

= Hi = \ {H'-' (x) + H'-' {x)) . (14) 

Hf = HZ = \ (iJ^=" {x) - H'=^ {x)) . (15) 

The pion valence parton distribution (ttVPD) is defined as Hy^ (x) = (x) + i?^ {—x) . Using charge conjugation, 
one can prove that the isoscalar and isovector part the quark distribution, H^^^'^ [x) , are respectively odd and even 
functions of x. Therefore 

H]}\{x)^\ J ^e^^- (tt' iP)\^iO) 7+i«r^^(An)|^MP)) 

+ IJ ^e-'^^ {7T^{P)\{p{0)j+^aT'^P{Xn)\7T^{P)) (16) 
In terms of the standard q and q distributions we have 

Hlf, {x) = ia i e,3, [qv (x) 9 (x) + qy i^x) 9 i-x)] (17) 

with qy (x) = q (x) — q{x) . 

Soft pion theorems have been used to relate the double pion distribution amplitude to the TrDAfsH- Using the 
same arguments, it is easy to connect the matrix elements of the preceeding definitions, 

]iTn{Tr'{P)\4,f{Q)^+T^^:f{Xn)\T:= {P))=ieai^ lim (0| ^/ (0) 7+75tV/ (An) Itt^P)) • (18) 



and therefore 



{x) = \{a)i e,3, {x) 9 (x) + 0^ {-x) 9 {-x)] . (19) 



2 

Substituiting Eq. in Eq. we obtain 

(f>^ix) = qvix), (20) 

a relation which is important for our scheme because it provides a universal determination of the hadronic scale Qo- 
The previous result, Eq. (PD|) . is confirmed by chiral models for the pion, like the NJL model or the spectral quark 
model. In fact, these local chiral quark models give, in the chiral limit, (j>.n{x) — q{x) — 1. Using the results obtained 
in the NJL model for the nDA (see [l3|) and the ttPD (see [s^) we obtain 

<PAx) = l + 0{ml) , qy(:,) = 9^ + o{ml)=l + 0{ml) (21) 
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where in the last relation ruq is the constituent quark mass, g^rgg is the pion quark coupling constant and the last step 
follows from the Goldberger-Treiman relation. This result for qvix) arises naturally from the physical meaning of 
the ttVPD. The flatness of the distribution implies that the probability of finding a quark of any relative momentum 
fraction x is the same. This is so because in the chiral limit quarks and pions are massless and therefore there is no 
mass scale to which the momentum distribution can attach. The appearance of flat distributions is a fundamental 
consequence of soft pion theorems and the chiral limit, not a peculiarity of some models. 

Since the origin of the flat distribution is related to soft pion theorems and the chiral limit, this behavior will change 
as we move towards higher energy. Starting from Qq, the momentum scale in which the flat distributions are a good 
approximation, 4'tt{x,Qq) = qv{x,Qo) — 1, we apply the evolution equations in order to obtain these distributions 
at higher Q values. The nVPD, qy{x,Q) evolves through the DGLAP equation, and it concentrates close to the 
a; ~ region for very high values of Q. The ttDA evolves through the ERBL equation and, therefore, for high Q 
values we obtain the asymptotic result (I)t^{x,Q) ~ Qx{l — x). In conclusion, QCD evolution obscures the relation 
between these two distributions, which become distinct from one another. It follows from the above discussion that 
it is unreasonable to assume that at relatively low energies <1)tt{x) is close to its asymptotic value, whereas the nVPD 
is far from it. 

The evolution of the flat ttVPD, q{x,Qa) — 1, has been previously studied in refs. [13, [2l| - [23j obtaining a very 
good description of the experimental data. We sketch here the procedure to fix the value of the matching point, Qo- 
We know that the momentum fraction of each valence quark at Q = 2 GeV is 0.23 [s^l- We fit the initial point of the 
evolution, Qqi imposing that the evolved irVPD reproduces this momentum fraction at Q = 2 GeV. Obviously, the 
resulting Qq value will depend on the value of Aqcd- Varying Aqcd form 0.174 GeV to 0.326 GeV we obtain that Qq 
changes from 0.290 GeV to 0.470 GeV, in leading order (LO) evolution. A change in Aqcd implies a change in the 
initial point Qq but, once we have fitted Qq through the momentum fraction carried by the quarks at Q = 2 GeV, the 
evolved parton distribution is independent of the value of Aqcd and the quality in the description of the experimental 
data is preserved. 

At this point we must return to Eq. (fTO| regarding the stability of AI. We obtain that M changes only a few MeV 
when we vary Aqcd form 0.174 GeV to 0.326 GeV. The largest incertitude comes from the experimental error in a, 
obtaining M — 0.466±0.006 GeV . We thus obtain a very stable value for M, showing the consistency of the procedure. 
For the numerical analysis we have chosen Aqcd = 0.226 GeV, Qq = 0.355 GeV, a — 0.32 and M = 0.466 GeV . 

IV. DISCUSSION. 

Summarizing the ideas developed previously, a model independent scheme of evaluation for the ttTFF at high Q 
values arises under the following assumptions: i) we consider, on the basis of chiral symmetry and soft pion theorems, 
that at some point Qq the ttDA and the nVPD coincide and have flat behavior (j)-„-{x,QQ) — gy(a;, Qo) = 1; **) we 
fix Qo applying QCD evolution to the ttVPD from Qq to Q = 2 GeV and impose that each valence quark carries a 
fraction of the total momentum of the pion of 0.23; Hi) for Q < Qq we assume the experimental parametrization of 
F{Q^) given in ([H]); and iv) for Q > Qq the ttTFF is given by the Eq. where we fix the M parameter through 
the continuity condition Eq. ([TUl) . 

In Figs. [2] and [3] we show the ttTFF calculated using the evolved irDA as explained above and compare with the 
data of refs. \X,M,i2Z]- 

In ref. [2l| for the pion, and in refs. [l^ for the nucleon, it is shown that including gluons and sea quarks at 
low energy the preferred starting point for evolution is Qo — 1 GeV. To move Qo up to 1 GeV implies to change the 
value of M to 0.500 GeV. Dashed curves in Fig. [2] and [3] reproduces the ttTFF calculated with these parameters. 
The experimental incertitude in a produces in this case a larger indetermination in the mass, obtaining M = 0.500 ± 
0.040 GeV . Nevertheless, this error does not affect to the predictive power, as shown in Fig. [51 

The calculation shows several features which we now discuss. The first effect of this model independent calculation, 
in comparison with the triangle result, is to cure the asymptotic behavior of the form factor. Model calculations 
overshoot the experimental data for ~ 10 — 30 GeV^ as is clear from ref. pjj and it is confirmed by our calculation 
of the triangle diagram (see appendix), where we obtain Q'^F{Q^) ~ 0.54 GeV for Q^ ^ 30GeV^. Moreover, the 

form factor in this model calculation shows a behavior proportional to [log (Q^/to^)] for high Q^. Now, with the 
procedure defined here, we recover the right asymptotic result Q^F (Q^) V^/rr- What is under discussion at present 
is not this value, but for which value of Q^ the asymptotic behavior is attained and the behaviour of the ttTFF in 
the intermediate energy region. 

In Fig. [3] we observe that the calculation with the mass cut-off introduced in ref. adjusts well to the data up 
to Q2 _ l5GeV^ when we use as matching point at Qo — 0.355 GeV. However, at that momentum most of the new 
data bend upward while our calculation remains practically fiat (see Fig [2]) . The slope of the ttTFF in the region 
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FIG. 2: Calculation of the transition form factor via the ttDA with M — 0.466 GeV, a = 0.032 and defining the matching 
point at Qo = 0.355 GeV (full-line) compared with the experimental data. The dashed-line corresponds to the same calculation 
defining the matching point at Qo = 1. GeV and with M = 0.500 GeV. The dotted-line corresponds to the asymptotic limit. 
The gray region gives the indeterminacy of the results due to the error in a, which is only apreciable for the Qo = 1 GeV case. 
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around 20 GeV^ is not well reproduced. If we assume the higher matching point Qo = 1 GeV, the calculation still fits 
relatively well the low energy data up to 5GeV^, but, contrary to the other, fails in the intermediate region, as we 
see in Fig. [3l In the region around 20 GeV^, this calculation goes through most of the experimental points although 
the slope of the curve is smaller than the one indicated by the data. The indeterminacy in the predictions related to 
the experimental error in a is only apreciable when we use Qq — 1 GeV as matching point. It is represented by the 
gray region in Fig [21 corresponding the upper limit of the region to a = 0.028 and the lower limit of to a = 0.036. In 
both calculations, the ttTFF tends to the asymptotic limit at very high energy, but in the studied region it is growing 
very slowly {Qq = 1 GeV case) or it is practically constant {Qq = 0.355 GeV case). The value of the mass parameter, 
around 0.450— 0.500 GeV, is reasonable, since it receives contributions not only from de constituent quark mass, but 
also from the mean value of the transverse momentum. 



V. POSSIBLE TWIST THREE EFFECTS IN THE PION TRANSITION FORM FACTOR. 

Using Qq = 0.355 GeV as matching point, the experimental results are well reproduced up to 15GeV^, but it is 
difficult to accommodate the higher momentum results. QCD is assymptotically free, and the perturbative approach 
leading to Eq. ([l} must be right for high enough momentum transfer. To correct the high momentum region without 
touching the low momentum region is not easy, and might require new effects, as the one shown in ref. 6]. 

On the other hand, using Qo — 1 GeV as matching point, the high momentum results are well reproduced, but with 
a too small slope, and the lower momentum data are over estimated. Here we are in a more confortable situation from 
the point of view of QCD. In order to cure the divergences in the integrand of Eq.([l]) we have introduced a cutoff mass 
M in the definition of the ttTFF in terms of the ttDA, recall Eq. ([6]). This procedure incorporates effects of twist 
three operators at low energy. The set of twist three distribution amplitudes for pions are (0| J (0) i 75 r* u (z) [tt^^ , 
(0| d (0) cTafi 75 T^u (z) [tt-') and (0| d (0) ctq^ 75 r* Gfj,u (y) u (z) [tt-') . Some of the effects of these operators are included 
in the modification proposed by Radyushkin, but there can be additional effects, like the presence of gluons, which 
are not incorporated in his procedure. These additional effects can be introduced by adding to the lowest order 
calculation a term proportional to , 

Q'F{Q') = ^f -^M^, Q') + §- (22) 
^ Jo X + -qt V 

The effect of C3 decreases rapidly with Q^, and therefore we consider C3 constant in Q^. Even though this constant 
is connected to twist three operators, its QCD evolution is unimportant here. 

Following the same ideas developped in Section [ill we change the matching condition, Eq ([TU)) in order to introduce 
the effect of the new term: 



V2f., Ql + M^ , C3 f (0) Ql 

With Qo = 1 GeV, this equation allows to determine M, once we have fixed the value of C3 . 

We consider three scenarios, which correspond to a contribution from the twist three term to the form factor at 
Qo = 1 GeV of 10% (C3 = 0.99 lO'^ GeV=^), 20% (C3 = 1.98 10"^ GeV^) and 30% (C3 = 2.98 lO'^ GeV^). If we fix^ 
C3 = 1.98 • 10~^ GeV^, which implies that this term is responsible for a 20% of the value of the nTFF at Qo = 1 GeV, 
the experimental error in a produces a variation of M = 0.620 ± 0.050 GeV. To check the stability of the M value 
in relation to C3, we fix a = 0.032 and vary C3 from C3 = 0.99 10"^ GeV^ to C3 = 2.98 10"^ GeV^ . This produces 
a variation of M given by M = 0.620 ± 0.070 GeV. The indeterminacy in M is in all cases of the order of 10%, 
which is of the same order of the relative error in the experimental determination of a. The value of M is sufficiently 
stable to confirm the consistency of the procedure. For the numerical analysis we have chosen Aqcd = 0.226 GeV, 
Qo — I - GeV, a — 0.032 and C3 — 1.98 10~^ GeV^ which corresponds to a 20% contribution of twist 3 at Qo- This set 
of parameters implies M = 0.620 GeV. 



^ A dimensionless definition of the parameter C3 could be C3 = 03 ym'^g/aj /■„-, with a the experimental parameter given in Eq. SEl . 
Therefore, C3 = 1.98 10~^ GeV"^ corresponds to a3 = 0.27. In the Appendix we give the value for C3 obatined from the triangle diagram. 
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FIG. 4: Calculation of the transition form factor via the nDA with M = 0.620 GeV, a = 0.032, C3 = 1.98 10"^ GeV^ and 
defining the matching point at Qo = 1 GeV (full-line) compared with the experimental data. The gray region gives the 
indeterminacy of the results due to the indeterminacy on C3. 

In Fig. |3]and[n]we show the ttTFF calculated using the evolved nDA compared with the data of [ll[H,[23|. In Fig 
[5] we have used the experimental parametrization ([8]) for < 1 GeV^ (short dashed line). For > 1 GeV^, the full 
line corresponds our calculation using Eq. (I22p together with the evolved ttDA. The gray region corresponds to the 
effect due to the different contributions of the C3 term: the upper limit corresponds to C3 = 0.99 IQ-^GeV^ (10% of 
twist 3 contribution at = 1 GeV); the lower limit corresponds to C3 = 2.98 10~^ GeV'^ (30% of twist 3 contribution 
at = 1 GeV); the dashed curve in Fig. [5] corresponds to the contribution coming from the C3 term for our central 
value, C = 1.98 • 10~^ GeV"^ . We observe that its contribution is absolutely negligible for > 5 GeV^ . Nevertheless, 
its influence in the matching point is important for a good description of the low data. Looking at the whole of 
the Figs. |3]and[5l we observe that we have a good determination of the experimental data in the whole experimental 
region. 

VI. CONCLUSIONS 

In the previous sections we have developed a formalism to describe the ttTFF on all experimentally accessible 
range, and hopefully beyond. The formalism is based on a two energy scale description. The formulation in the 
low energy scale is non perturbative, while that of the high energy scale is fundamentally a QCD based perturbative 
formulation. The two descriptions are matched at an energy scale Qq hereafter called hadronic scale. This scheme has 
been applied to described the parton and generalized parton distributions with notable success [l8l - l20l . [ssj . In order to 
understand the finesse of the BaBaR data we have forced us to introduce additional ingredients to this scheme which 
we next summarize. 

We have taken for the low energy description of the ttTFF a parameterized description of the data. This has 
been done to avoid model dependence in a region (Q ~ 1 GeV), where models might start to break down. The high 
energy description incorporates the following important physical ingredients: a mass cut-off in the definition of the 
ttTFF from the ttDA, M, which, interpreted from the point of view of constituent models, takes into account the 
constituent mass, transverse momentum effects and also higher twist effects; ii) a flat irDA at low energies in line with 
chiral models pH - flsl . [ist and recent proposals 0,0] ; iii) eventually an additional twist three term into the definition 
of the ttTFF in the high energy description parameterized by a unique constant C3; iv) the two descriptions have to 
match at an energy scale Qq; this scale is universal and should be the same for all observables. Though our scheme 
is not based in the use a hadron model, we have presented in the Appendix a simple model calculation to illustrate 
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FIG. 5: Same as in Fig |4] for the region under 17GeV'^ . The dashed-hne corresponds to the C3 contribution. The short 
dashed-hne is the experimental parametrization ((Sj for < 1 GeV^ . 

the contents of the procedure. 

We have seen that, as a consequence of soft pion theorems, the ttDA and ttVPD are equal at low energies up to the 
matching scale Qo- Moreover, chiral models indicate that in the chiral limit both distributions are flat, more precisely, 
07r (x, Qa) = qv {x, Qo) = 1- The flatness of the ttVPD, which implies that the probability of flnding a quark of any 
relative momentum fraction x is the same, is justified by the fact that in the chiral limit quarks and pions are massless 
and therefore there is no mass scale to which the momentum distribution can attach. The different tipe of evolution 
followed by these two distributions hides the relation between them. Nevertheles, the equality between the ttDA 
and the ttVPD at Qo scale allows us to determine the matching scale using the parton distributions. From previous 
studies [13, m, [2^ we can fix this scale at Qo = 0.355 GeV. Nevertheless, more elaborated studies (l9l - l2T| indicates 
that this matching scale cannot be one in which the constituent quarks are purely valence but require additional sea 
and gluon components. Theses studies push the matching scale at Qo ~ IGeV. 

Using Qo — 0.355 GeV as matching point, the experimental results are well reproduce up to 15GeV^, but it is 
difficult to accommodate the higher momentum results. If we adopt this philosophy, the BaBar data, if confirmed, 
might imply some new type of contributions @. Using Qo = 1 GeV as matching point, the higher momentum 
results are well reproduced, but with a slope which is too small, and the lower momentum data are overestimated. 
Nevertheless this last situation can be corrected in a natural way if we assume some additional contribution from 
twist three. 

The calculations shown prove that the BaBar results can be accommodated in our scheme, which only uses standard 
QCD ingredients and low energy data. It must be emphasized that in order to have a good description higher twist 
effects are important as the modification from Eq. ([T]) to Eq. signals. It must be also noted that the matching 
scale is as high as 1 GeV, a feature which was also the case in the description of parton distributions when precision 
was to be attained. With these ingredients our calculation shows an excellent agreement with the data and points 
out that maybe the average value of the highest energy data point is too large, a conclusion reached by other analyses 



The relative high value of M — 0.620 GeV can be understood on the basis that it includes the constituent quark 
mass, the mean value of the transverse quark momentum and other higher twist contributions. The C3 term is 
relatively small. Its effect is to reduce the value of the contribution to the ttTFF of twist two only for low Q^ < 5 
GeV^. Our results are very stable with respect to variations of these parameters. 

Let us conclude by stressing that we have developed a formalism to describe the ttTFF based on the philosophy 
that two descriptions, a perturbative and a non perturbative can be matched in a physically acceptable manner at 
a certain scale Qo. The idea of the formalism is that one can use models or effective theories to describe the non 
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perturbative sector while QCD to describe the perturbative one. In here we have preferred to use data for the low 
energy sector to avoid model dependence, but in other observables, and specially in predictions, models can be used 
for estimates. Moreover, in order to describe the data we have seen that a matching scale of 1 GeV was favored over 
the one of 0.355 GeV that we would obtain if only valence constituents were assumed to describe low energies. This 
implies that the description of high precision data requires sophisticated effective theories (models) at low energies. 
We have also discovered that higher twist effects (parametrized in our case by M and C3) are small but crucial in 
order to attain good descriptions. Finally we can assert that our description of the data is excellent all over although 
points out to a too large average value of the highest Q points in the data. 
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Appendix A: The triangle diagram. 



The easiest model to evaluate the ttTFF is to use the triangle diagram defined in Fig. [TJ There are two diagrams, 
the one depicted an another with the photons exchanged. The amplitude for the process can be defined as 



In a straightforward calculation"^ 



where 



2 „'2\ 



r = 4^a£^e"'(7- q"' e^.p. F {q,q') 



F{q^,q''^) = 8mqg^qql{m,q^,q'^), 



I{m,q ,q 



(27r) (^P-m^^+ie) {q-tf -ml+ie {q' + if - + i 



(Al) 
(A2) 

(A3) 



qfj, and q'^ are the photon momenta, which, in this expression, can be real or virtual photons. Here, mq is the 
constituent quark mass, g^^qq the pion-quark coupling constant, q and g' the photons momenta. The small virtuality 
of one of the photons implies that q'^ — 0, and we define F{Q^) ~ F{q^,0) with — —q^. In the chiral limit, = 0, 
we can use the Goldberger-Treiman relation, g-^qq/mq = \/2//7r, obtaining 



\/2m 



2 7r2/. 



dx- 



1 



log 



1+M + VTT2I 



1 + u - VT+2u 



(A4) 



where u — Irn^ J {Q"^ x^) . 



The behavior of the form factor for small values follows from the previous equation. 



1 - 



12 m2 



(A5) 



m^/\/T2^~ 0.220 GeV. 
41 can be developped in the form 



Therefore, we have F(0) = V2/ {'in^U) = 0.273 GeV~\ and ,from Eq. ©, m 
We can also look for the asymptotic behavior. The denominator of Eq. 

xVl + 2 M ^ X + 2mq/{Q'^x). Comparing this expression with Eq. ©, we have M — ^2m'^/x ~ 2mq = 0.44 GeV 

Performing the integral in Eq. (jA4p we have 



Q2f(Q2) 



>oo 2V27T^fr 



1 

■log — 



log — 



n2 



(A6) 



^ A detailed calculation of this process is shown in ref. [35| . We need to replace there by Nc (ej — e^) = e^, since in our case we have 
quarks in the loop . 
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Therefore, the model calculation gives a wrong asymptotic behavior. We have Q'^F{Q'^) ~ 0.54 GeV for ~ 30 GeV^ 
which is a very large result, and the form factor grows as (log Q^/mg) for large values. Comparing this expression 
with Eq. ((22|) . we obtain C3 = v^rn* log ((5o/™g) /""^/tt ~ 0.78 10~^ GeV^, which is consistent with the values for 
C3 used in the text. 
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